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ABSTRACT 

The emission mechanism responsible for the bulk of energy from radio to X-rays in low ionization emission line regions 
(LINERs) and Low Luminosity Active Galactic Nuclei (LLAGN) has been long debated. Based on UV to X-ray and 
radio to UV flux ratios, some argue that LINERs/LLAGN are a scaled-down version of their more luminous predecessors 
Seyfert galaxies. Others, based on the lack of X-ray short (hours) time-scale variability, the non detection of an iron 
line at 6.4 keV, and the faint UV emission compared to typical AGNs, suggest the truncation of the classical thin 
accretion disk in the inner regions of the AGN where a radiatively inefBcient accretion flow (RIAF) structure forms. 
We investigate the LINER-Seyfert connection by studying the unabsorbed LINER galaxy NGC 4278 that accretes at 
a low rate (Lbot/Edd ~7xl0~'') but exhibits a broad Ha line, and a point-like nucleus in radio, optical, UV and X-rays. 
We analyzed one XMM-Newton and seven Chandra X-ray observations of NGC 4278 spread over a three year period, 
allowing the study of the X-ray variability at different time-scales (hours, months, years). We also examined the radio 
to X-ray spectral energy distribution to constrain the accretion mode in the nucleus of NGC 4278. Long time-scale 
(months) variability is observed where the flux increased by a factor of ~3 on a time-scale of a few months and by 
a factor of 5 between the faintest and the brightest observation separated by ~3 years. During the XMM-Newton 
observation, where the highest flux level is detected, we found a 10% flux increase on a short time-scale of a few hours, 
while the light curves for the different Chandra observations do not show short time-scale (minutes to hours) variability. 
A combination of an absorbed power law {Nh ~ 10'^'' cm~^, F = 2.2j;Q2) plus a thermal component (kT~ 0.6 keV) were 
able to fit the Chandra spectra. The XMM-Newton spectra, where the highest X-ray flux is detected, are well fitted 
with an absorbed power-law with no need for a thermal component as the emission from the power-law component is 
dominant. The power-law photon index is ~ 2.1 and the hydrogen column density is of the order of 10^° cm~^. Neither 
a narrow nor a broad Fe Ka emission line at 6.4 keV are detected with a 22 eV and 118 eV upper limits derived on 
their equivalent widths. We derive optical fluxes from archival HST ACS observations and detected optical variability 
on time-scales of years. For the first time for this source, thanks to the optical/UV monitor on board XMM-Afewton, 
we obtained simultaneous UV and X-ray flux measurements. We constructed SEDs based on simultaneous or quasi 
simultaneous observations and compared them to LINER, radio-loud, and radio-quiet quasar SEDs. We find that at 
a low X-ray flux the NGC 4278 SED resembles that of typical LINER sources where the radio to X-ray emission can 
be considered as originating from a jet and/or RIAF, whereas at a high X-ray flux, NGC 4278 SED is more like a low 
luminosity Seyfert SED. Consequently, NGC 4278 could exhibit both LINER and Seyfert nuclear activity depending 
on the strength of its X-ray emission. 
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1. Introduction 

Low-ionization nu clear emission line regions (LINERs) were 
first identified by Heckman ( 1980 1 as a class of galaxies 



with optical spectra dominated by emission lines from low 
ionization species. The ionization mechanism is yet poorly 
known and could be explaine d, either in terms of starbursts 
Alonso-Herrero e t al.||2000[) or, more consistently, being 
ue to a low luminosity active gala ctic nuclei (LLAGN) 



( Terashima et al.|2000( ,Ho et al. 1993 1 . This latter idea rises 
from the detection of broad Ha components in an important 
fraction of LINE R sources (noted as LINER 1.9 (hereinafter 



LINERI) objects Ho et al. 1997), and/or a point-fike UV 



or X-ray source at the nucleus. If the luminosity scales with 
the accretion rate, the study of the LINER nucleus define 



a supreme pattern for probing low accretion rate physics 
around supermassive black holes (SMBHs). 

How similar are these LINERI sources to classical lumi- 
nous Seyfert and quasar galaxies? Based on observational 
properties, the weakness or absence of a big blue bump fea- 



CMalkan & Sargent||l982 Sanders et al. 


19891 


Koratkar & 


Blaes 1999), the lack of a broad Fe Ka emission line at 


6.4 keV (Terashima et al. 

j_l _ TT 1 1 M 1 11 S \W 


2002; Ptak ct al. 2004|), except 



and last but not least, the lack of short time-scale 
utes to hours) X-ray variability have been attributed to 
an intrinsic difference in LINERsI/LLAGN central engine 
as opposed to normal Seyfert galaxies. One proposed sce- 
nario is that accretion in LINERsI/LLAGN is radiatively 
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inefHcient, advection dominated compared to the typical 
geometrically thin optically thick accretion disks present in 
luminous AGN (Narayan. 2005) . 



However, more recently jMaoz et al. ( 2005 ) revealed UV 
variability in the nuclei of a sample of 17 LINERs ob- 
served with the HST. By combining these res ults with non- 
simultaneous radio and X-ray observations, Maoz (2007) 



demonstrated that the UV/X-ray luminosity ratios are 
similar to those of Seyfert 1 nuclei and pointed out that 
LLAGN may be a scaled-down version of Seyfert galaxies 
where a thin accretion disk exists. This idea is supported 



by Plan et al. (2010) who studied simultaneous UV to X- 
ray observations of four LINER nuclei observed with the 
XRT and the UVOT on-board the Swift telescope. They 
discovered short-time scale (half a day) X-ray variability in 
two of their objects and showed that the UV to X-ray flux 
ratios are consistent with those of more luminous AGN. 

The elliptical g alaxy NGC 4278 (distance of 16.7 Mpc; 
Tonry et aL||200lj scaled to Ha = 70 km s^^ Mpc^i) has 



been studied extensively at different wavelengths. At ra- 



been detected at 6 



Jones et al. 


1984 


1. 


Nagar 



sec scales in the form of twin jets emerging fr om a cen- 
tral co mpact component {Tb — 1.5 x 10^ K). Ho et al. 
" 19971 classified NGC 4278 as a type 1.9 LINER from the 
ehnite detection of a broad Ha line, with log F{Ha) = 
— 13.07 ergs cm~^ s~^, supporting the AGN nature of 
the nuclear engine. Moreover, an unresolved compact nu- 
clear source has be en detec ted in HST WF PC2 images 
( Capetti et al.|2000| . Finally, |Ho"et al.| ( |2001[ ), after study- 
ing a Chandra 1.4 ks snapshot taken in April 2000, gave 
NGC 4278 a class ( I) X-ray morphology, show ing a domi- 
nant nuclear source. Terashima & Wilson ( 2003 1 fit the 0.5- 
8 keV spectrum obtained during this same Chandra snap- 
shot with a power law modified by absorption and found 
a 2-10 keV corrected luminosity of 9.1 x 10'^^ ergs s~^. 



luminosity of 5.6 



Gonzalez-Martm et al. (2009) find a 0.5-10 keV corrected 



400 ks 



X iU"" erg s"" after studying a 
Chandra observation taken in March 2006. 

The mass of a black hole in the nucleu s of NGC 42 78 
derived from the M- cr relation ([Tr cmainc et al. 2002 ), is 



3.09z b0.54xl0^ Mq ( |Wang fc Zhang|2003f |^iaberge et al 
2005 1 . This leads to a LEdd » 3.9 x 10'^° ergs s~^; assuming 
tHaFXsdd = 1-25 X 10^^ (Mbh/Mq) ergs s'^ This implies 
a very low Lx/^Edd ratio of 2xl0~^. 

In this paper we report a timing and X-ray spectral 
study of the LINER galaxy NGC 4278 observed with XMM- 
Newton and Chandra. We describe in Sect. 2 the X-ray ob- 
servations and the data reduction. Timing and spectral X- 
ray results as wefl as optical (HST/ACS) and UV (XMM- 
iVewtoii/OM) results are presented in Sect. 3. In Sect. 4, 
we describe the construction of the NGC 4278 spectral en- 
ergy distribution. We discuss in sect. 5 our results in the 
context of the LINER-Seyfert connection. The main results 
are summarized in Sect. 6. 



taken over a two year interval starting on 2005 February 03 
and ended on 2007 April 20. The seven observations were 
obtained with the spectroscopic array (ACIS-S; Weisskopf] 
|et ar.|2002| where the nucleus was placed on the aim point 
of the ACIS-S3 back-illuminated chip. They were taken 
in Very Faint mode to increase their sensitivity. All of 
the observations are Chandra archival data obtained from 
chaser ^ . The log of the Chandra observations are listed in 
Table 1. 

All seven observations were reduced and analyzed in 
a systematic, homogeneous way using the CIAO soft- 
ware package version 4.2, Chandra Calibration Database, 
CALDB, version 4.2, an d the ACIS Extrac t (AE) software 
package version 3.175 2 (iBroos et al.||2010|). We started by 
using the level 1 event file produced by the Chandra X- 
ray Center (CXC) to suppress the position randomization 
applied by the CXC Standard Data Processing when cre- 
ating a level 2 event file. We also corrected for the effects 
of charge-transfer inefficiency on event energies and grades. 
We filtered for bad event grades (only ASCA grades 0, 2, 3, 
4 and 6 are accepted) and hot columns to take account of 
several potential issues such as cosmic rays and bad pixels. 
Good time intervals, supplied by the pipeline, were applied 
to the final products 

In order to get the X-ray position for the sources in the 
ACIS-S field of view, we used a wavelet transform detec- 
tion algorithm, the w avdetect program with in the CIAO 



data analysis system (Freeman et al. 2002). The output 



file of the previous procedure is input to the AE software 
The AE package is used to refine source positions, ex- 
tract source photons, construct local backgrounds, extract 
source, and background spectra (see below), compute redis- 
tribution matrix files (RMFs) and auxiliary response files 
(ARFs), by spawning the mkarf and mkacisrmf routines 
of CIAO, and perform spectral grouping and fitting. We fo- 
cus here on the X-ray emission from the nuclear source (for 
inf ormation on the X-ray s ource population of NGC 4278 
Brassington et al.||2009[ ). 



see 



Source events are extracted around the source centroid, 
inside a polygonal shape of the local PSF, generated by 
MARX'' at the energy 1.497 keV using the ae-make-psf 
tool implemented in AE. Background region was defined 
on the brightest of the seven observations (obs ID: 4741) 
following the AE procedure. This background region is an 
annular region centered on the source position where the 
inner radius delimit 1.1 x 99% encircled energy radius and 
the outer radius is set such that the background includes 
200 counts. This number is obtained from a special im- 
age where all events within the ~ 1.1 x 99% PSF circles 
of all the sources in the field were excluded {swiss cheese 
image). This background region is then generated to the 
other observations in order to have a consistent and coher- 
ent measurement of the diffuse/unresolved X-ray emission 
between the seven observations for a homogeneous spectral 
analysis. 



2. X-ray observations and data reduction 

2.1. Chandra observations 

NGC 4278 was observed seven times with the Chandra oh- 
ser vatory, with a fir st 1.4 ks snapshot taken in 2000 April 
20 (Ho et al. 2001). The following six observations were 



http: / /cda.harvard.edu/chaser/Descriptions 



The 
and 



User's 



Extract 
Guide 



softwaxe 
are available 



pack- 
at 



http:/ /www. astro. psu.edu/xray/acis/acis_analysis. html 

Description on how to create a new level 2 event file can be 



found at http://cxc.harvard.edu/ciao/threads/createL2/ 
http://space.mit.edu/ASC/MARX7 
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Table 1. Log of the X-ray observations. 



Satellite 


Instrument 


Start date 


Obs. ID 


Exposure time (ks) 


Chandra 


ACIS-S 


2000 April 20 


398 


1.4 


XMM-Newtoif 


EPIC 


2004 May 23 


205010101 


34/35.5/35.5 


Chandra 


ACIS-S 


2005 February 02 


4741 


37 


Chandra 


ACIS-S 


2006 March 16 


7077 


110 


Chandra 


ACIS-S 


2006 July 25 


7078 


51 


Chandra 


ACIS-S 


2006 October 24 


7079 


105 


Chandra 


ACIS-S 


2007 February 20 


7081 


110 


Chandra 


ACIS-S 


2007 April 20 


7080 


55 



Note: "Exposure time for pn/MOSl/MOS2 respectively. 



We find that all of the seven Chandra observations 
suffer from pile-up. X-ray data affected by pile-up can 
mimic spectral and timing variability. Therefore, pile-up 
was accounted for by excluding the core of the PSF (see 
Appendix A for more details). We decided not to use the 
first 1.4 ks Chandra snapshot observation of NGC 4278. 
Indeed, after removing the piled-core, not enough counts 
remained to build a useful spectrum. We use the tool dmex- 
tract, called by the AE software, to create a source spec- 
trum over the energy range 0.5-8 keV. We used the tool 
ae^group^spectrum implemented in AE to group the spec- 
tra. Channels between 0.5 and 8 keV are grouped to have 
a three sigma {3a) signal to noise ratio, which corresponds 
to a minimum of 20 counts per bin, to enable the use of 
statistics in the spectral analysis. 

2.2. XMM-Newton observation 

The log of the XMM-Newton observation is listed in 
Table 1. NGC 4278 was obser ved by XMM-N ewton on 
2004 May 23. The EP IC-pn ( jStriider et al.|2001[ ) and MOS 
(Turner et al. 2001) cameras were operated in Imaging, 
Prime Full Frame mode using the thin filter and medium 
filter respectively (exposure time about 34 and 35.5 ks, re- 
spectively) . The Reflection Grating Spectrometers show few 
counts therefore they were not included in our analysis. 
The optical/UV monitor (OM) instrument (Mason et al. 
2001 1 was operated in imaging mode with the UVWl and 



the UVM2 filters (effective A of 2910 and 2310 A, respec- 
tively). All data products were obtained from the XMM- 
Newton Science Archive (XSA)^ and reduced using the 
Science Analysis System (SAS) version 9.0. The EPIC data 
were not affected by pile-up, and no severe intervals of en- 
hanced solar activity were present. Data are selected using 
event patterns 0-4 and 0-12 for pn and MOS, respectively, 
during only good X-ray events ( "FLAG=0" ) . 

We extract source events from a circle of 10" radius 
centered on the source. Background events are extracted 
from a source-free circle of the same radius on the same 
CCD as the source. We generated response matrix files us- 
ing the SAS task rmfgen, while ancillary response files were 
generated using the SAS task arfgen. The EPIC spectra 
were created in the energy range 0.5-8 keV to enable flux 



and model-parameter comparison with Chandra spectra^. 
They were grouped using the FTOOLS grppha tool to have 
a minimum of 20 counts per bin to allow the use of the 
statistic. The SAS task omichain was used to process the 
OM data. 



3. Results of the nucleus study 

3.1. Lightcurves and hardness ratios 

Light curves for the six Chandra observations were ex- 
tracted in the energy range 0.5-10 keV from the pile-up 
free wings of the PSF. To enable the comparison between 
the different observations, the net count rate of each is then 
corrected for the excluded fraction of the PSF. The hard- 
ness ratio is defined as, HR = {H — S)/{H + S), where 
S is the count rate in the soft 0.5-2 keV band and H is 
the count rate in the hard 2-10 keV band. Fig. [l] shows the 
0.5-10 keV light curves and the hardness-ratios of the differ- 
ent observations. No clear-cut short timescale, minutes to 
hours, variability is detected in any of the observations. We 
conducted a Kolmogorov-Smirnov, K-S, test to check more 
accurately any variability in the source count rate within 
each observation. We find a K-S test probability > 17% 
in the six observations that the nuclear emission originates 
from a constant source. Therefore, NGC 4278 is considered 
as constant on short (minutes to hours) timescale. On the 
other hand, long timescale (months) variability is observed 
between the different observations. The hardness ratio in all 
of the six observations indicate that the source is in a soft 
state being constant during each observation and between 
the observations as well implying similar spectral charac- 
teristics. 

We extracted an EPIC-pn light curve from a 10" circle 
centered on the source. We used the tool epiclccorr'' to cor- 
rect for good time intervals (GTIs), dead time, and for back- 
ground subtraction. On account of the low angular resolu- 
tion of XMM-JVewton compared to Chandra, events from 
the 10" circle are contaminated by emission from point- 
like sources and diffuse emission (Fig. [2|. Nevertheless, we 
show in Appendix B that emission from the central source 



http: / /xmm.esac.esa.int/xsa/index.shtml 



^ In fact, the contribution from the surrounding medium 
around NGC 4278 that we calculated using the Chandra im- 
ages is unknown below 0.5 keV as we extracted the spectra in 
the 0.5-8 keV range, see sect. 3.3.2 for more details. 

http://xmm.esa.int/sas/current/doc/epiclccorr/epiclccorr.htnil 
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Fig. 1. Light curves of the six Chandra observations along 
with their hardness ratios defined as HR — {H ~ S)/{H + 
S), where S is the count rate in the soft 0.5-2 keV band 
and H is the count rate in the hard 2-10 keV band. Light 
curves were binned to have 50 groups with a signal to noise 
ratio greater than 10. The dashed lines show the averages 
of the count rates and hardness ratios. 




Fig. 2. Left panel: the XMM-Newton/MOSl imago of 
NGC 4278 with a 10" circle centered on the nucleus. Right 
panel: the Chandra image (obs ID: 7077) with a 10"-radius 
circle centered on the nucleus. The lower angular resolu- 
tion of XMM-Newton has blended the central source, X- 
ray point-like sources, and diffuse emission, detected with 
Chandra, in one point-like source. 



1.8^ 




-0.7 u 1 1 1 1 1 1 1 1 L 

0123456789 
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Fig. 3. Upper panel: EPIC-pn light curve of NGC 4278 
during the XMM-Newton observation with a time bin size 
of 1.5 ks. Lower panel: hardness ratio. The dashed lines 
show the average on the count rate and hardness ratio. 
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Fig. 4. The residuals in terms of sigma of an absorbed 
power-law fit to the different Chandra observations. Some 
features are clearly present in the soft 0.5-2 keV band sug- 
gesting emission from diffuse hot gas. 



with an error < 15% and ^20% for the Chandra and 
XMM-Newton observations simultaneously. This result is 
discussed in sect. 5.2. 



is dominant and all the surrounding sources only contribute 
in 2% of the total X-ray fiux. Figure [S] shows the pn light 
curve and the associated hardness ratio. According to a 
K-S test, NGC 4278 shows variability during the XMM- 
Newton observation with a K-S probability of 0.8% that 
the emission comes from a constant source. This variability 
is clear at the beginning of the observation where we detect 
a 10% increase over a period of ~1 hour. 

Finally, we calc ulated the normalized "exc ess variance" 



3.2. Spectral analysis 



The spectral analysis was performed using XSPEC ( Arnaud 



1996) version 12.5.1. The updated photo-electric cross sec- 
tions and the revised solar abundances of I Wilms et al.l 
(2000) are used throughout to account for absorption by 



neut ral gas. An absorbed galactic column density derived 



( 1999 1 for the seven light curves, m order to check more 



from|Kalberla et al.|(|2005l) (obtained with the W3NH tool^) 



a'^rr,^ d efined as in Nandra et al. (1997) and Turner et al. and fixed to 2.07 x 10^'^ cm~^ was applied to all the models 

Spectral uncertainties are given using Ax^ of 2.71, corre- 
sponding to 90% confidence for one interesting parameter, 
and to 95% confidence for upper limits. 



accurately any possible variability during each observa- 
tion. We find that cr™„ < 0.01 in all of the observations 
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Table 2. Best fit parameters of a mekal plus an absorbed power-law models for the Chandra data and an absorbed 
power-law model for the XMM-Newton spectra. 



Obs. ID 


Nh" 


r 


PI norm. 


kT" 


EM"' 


Flux'' 


CoiT. flux' 






lO'O cm-^ 




10"^ Photons keV^' cm" 


^ kcV 


10"' cm-=' 


10"" erg 


10-" erg s-i 


% 


4741 


<6.78 


2.13 [2.00-2.28] 


4.28 [3.92-4.71] 


0.62 [0.58-0.66] 


2.64 [2.27-2.99] 


18.1 [17.1-18.8] 


18.5 [18.2-20.1] 


94 


7077 




2.26 [2.16-2.39] 


1.82 [1.68-1.99] 






7.8 [7.5-7.9] 


8.21 [8.0-9.0[ 


83 


7078 




2.34 [2.22-2.47] 


4.21 [3.90-4.61] 






16.0 [14.9-16.2] 


16.5 [16.1-18.5] 


93 


7079 




2.38 [2.28-2.50] 


3.76 [3.54-4.07] 






14.1 [13.3-14.3] 


15.05 [14.7-16.4] 


90 


7080 




2.02 [1.84-2.20] 


1.10 [0.96-1.26] 






5.9 [5.5-6.1] 


6.0 [5.7-6.8] 


80 


7081 




2.12 [2.00-2.25] 


1.25 [1.14-1.29] 






6.2 [5.9-6.4] 


6.4 [6.2-7.1] 


81 








Reduced 


= 0.93 for 310 d.o.f. 










0205010101 


3.82 [3.05-4.59] 


2.05 [2.03-2.07] 


8.07 [7.90-8.24] 






33.8 [32.8-33.2] 


34.5 [34.2-34.8[ 


100 



Reduced = l-Ol for 487 d.o.f. 



Note: "Parameter linked for the different Chandra observations. 

''The emission measure (EM) of the mekal model, E]VI= J nenndV. 
'^Fluxes in the 0.5-8 keV energy range. 

''The power-law component fraction to the total corrected flux. 
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Next, we fit the six spectra simultaneously in order 
to identify the best physical model suitable for the six 
Chandra spectra and to better constrain the model param- 
eters. For this purpose, we tried three different models: 

f . Two power-laws affected by absorption, where the first 
power-law represents the soft 0.5-2 keV emission from 
the AGN and the second power-law represents the 2- 
10 keV compton up-scattered emission by a corona. 

2. A power-law affected by intrinsic absorption and a ther- 
mal component. The power-law represents the direct 
continuum emission from the AGN whereas the thermal 
component takes into account possible emission from 
unresolved X-ray binaries and/or supernovae remnants. 

3. A power-law and a thermal component affected by 
absorption and an ad ditional power-law component 
( Guainazzi et al.||2005 ). This model is similar to model 
(2) but includes the possibility of Compton-scattered 
X-ray emission. 

IVIodel (1) gives an acceptable fit with a reduced 
of 1.38 for 311 d.o.f. However, residuals around 1 keV per- 
sisted in the different spectra. We then tried model (2) leav- 
ing the intrinsic absorption component, the power-law, and 
the mekal parameters free for the six spectra. The tem- 
perature and the emission measure (EM) of the thermal 
component, as well as the hydrogen column density, were 
consistent within the error bars during the different obser- 
vations and therefore linked between each other. The fit 
gives a reduced of 0.93 for 310 d.o.f. The mean average 
of the photon index, F, is 2.2^q 2. The temperature of the 
thermal component is 0.62 ±0.04 keV. The intrinsic hydro- 
gen column density affecting the power-law has an upper 
limit of 6.78 x 10^° cm~^. According to an F-test the sta- 
tistical improvement in the fit between model (1) and (2) 
is highly significant. Finally, we tried model (3) used to fit 
the spectra of Compton-thin Seyfert 2 galaxies. Similar to 
model (2), the parameters of the thermal component and 
the hydrogen column density were linked for the different 
spectra. We find a reduced of 0.93 for 309 d.o.f. The 

absorption hydrogen column density has an upper limit of 

http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 8 x 10^" cm~^. The average value of the photon index, F, 



Fig. 5. Left panel: data and best -fit simultaneous model 
of an absorbed power-law plus mekal to the six Chandra 
observations. Right panels: the residuals of the fit in terms 
of sigma. The features in the soft band were accounted for 
by the thermal component. See the electronic edition of the 
Journal for a colour version of this figure. 



3.2.1. Chandra spectral analysis 

To obtain a first guess about the models that would be best 
suitable to the data, we first fit the six pile-up corrected 
spectra separately. An absorbed power-law model gives an 
acceptable fit in the six cases with a power-law photon in- 
dex, F, between 2.3 and 2.8. The intrinsic hydrogen column 
density has an upper limit of 1.2 x 10^^ cm~^. In five of the 
six observations, the soft 0.5-2 keV band shows residuals 
around ~1 keV suggesting th e possible presenc e of emis- 
sion from hot diffuse gas (e.g., Flohic et al.||2006 ). Figure [4] 
shows the residuals from the single absorbed power-law fit. 
We note that the only observation that did not show any 
evidence for residuals around ^1 keV is the one with the 
least net source counts (Obs. ID: 4741). The additio n of a 
thermal component {mekal model; Mewe et al. 11985 1 takes 
account of these residuals in the soft band and improves 
the fit quality in the six observations. 
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is 2.3 ± 0.2. The temperature of the thermal component 
is 0.63^0 05 keV. The addition of a second power-law did 
not improve the quality of the fit. The parameters of the 
thermal model and the power-law model did not change be- 
tween model (2) and model (3) while the hydrogen column 
density is better constrained using model (2) . Moreover, ac- 
cording to an F-test, the improvement of the fit by adding 
a power-law component is not significant. 

Therefore, we consider model (2), an absorbed power- 
law and a thermal mekal component, as best suitable for the 
six Chandra spectra and we show in Fig.[5]the data and the 
best fit model with the residuals in the right panels. Table 2 
gives the parameters of the best fit model, the uncorrected 
and the corrected 0.5-8 keV flux, as well as the power-law 
fraction to the total corrected flux. 



3.2.2. XMM-Newton spectral analysis 

The Chandra image of NGC 4278 (obs. ID: 7077) is shown 
in the right panel of Fig. [2] The central bright source, the 
NGC 4278 nucleus, is surrounded by a 10"-radius circle. 
The left panel is the XMM- JVewton/MOSl image with a 
circle of the same size centered on the source. The lower 
angular resolution of XMM-Newton compared to Chandra 
blended the resolved and unresolved point sources, as well 
as the diffuse emission, in one point source. In order to ex- 
clude the contribution from the contaminating sources, the 
six spectra of the two brightest objects, marked source 1 
and source 2 in the Chandra image, obtained in the en- 
ergy range 0.5-8 keV from each observation were fit si- 
multaneously with an absorbed power-law, assuming that 
these sources are L ow Mass X-ray Binaries (LMXBs) ( Kim 
& Fabbiano 2004). We then excluded the two sources and 
the central object and fit the rest, obtained in the energy 
range 0.5-8 keV, with a combination of a thermal compo- 
nent and a power-law assuming contribution from diffuse 
emission and unresolved point sources. The point sources 
and the diffuse emission do not show significant variabil- 
ity between the six Chandra observations and therefore we 
assumed them to be constant (detailed results on the two 
point-like sources and the diffuse emission are presented in 
Appendix B). When fitting the XMM-Newton spectra ob- 
tained from the 10" circle around the point-like source, we 
included the contribution from all of the above (Source 1, 
source 2, and the diffuse emission that contribute only to 
2% of the total flux, see appendix B for more details) , the 
rest is then attributed to the nucleus. 

We started the analysis by fitting the pn, MOSl, and 
M0S2 spectra simultaneously with an absorbed power-law 
model leaving the normalization free to take account for 
slight cross calibration uncertainties between the three de- 
tectors. The fit is noticeably acceptable with a reduced 
of 1.01 for 487 d.o.f. The photon index, F, and the normal- 
ization of the power-law are 2.05 ±0.02 and (7.73 ±0.17) x 
10^"* photons keV^^ cm^^ s~^ respectively. The intrinsic 
hydrogen column density is 3.82 ± 0.77 x 10^° cm~^. We 
note that this value is better constrained than the value 
derived from the Chandra fit. Then, we tried to fit the 
XMM-Newton spectra with the same models used on the 
Chandra data. The combination of two absorbed power- 
laws, model (1), did not give a better fit than one sin- 
gle absorbed power-law with a reduced of 1.02 for 484 
d.o.f. The addition of an unabsorbed thermal component to 
the absorbed power-law, model(2), gives an acceptable fit 



with a reduced oi 0.99 for 482 d.o.f. An F-test indicates 
that the probability for such an improvement (compared 
to the fit with a single absorbed power law) to occur by 
chance is 0.25. Moreover, the emission from the power-law 
component dominates the total nuclear emission to a 99% 
level even at low energies. Therefore, the thermal compo- 
nent is not required by the data and was not kept in the 
model. Model (3), where an additional absorbed power-law 
is added to model (2), did not ameliorate the quality of the 
fit with a reduced of 1.01 for 479 d.o.f. Therefore, the 
absorbed power-law model is favored for the XMM-Newton 
spectra. We give in Table 2 the best fit parameters. Fig. |6] 
shows the data, folded model, and the fit residuals. 

No evidence is found for an Fe Ka emission line at 
6.4 keV, neither in the pn spectrum nor in the Chandra 
spectra. The upper limit derived on its equivalent width 
(EW) is measured from the XMM-Newton observation as 
the spectral resolution of the EPIC-pn is higher than that of 
the ACIS-S Chandra instrument. By including a Gaussian 
profile with a width fixed to and 0.3 keV (for a nar- 
row and a broad line respectively) at 6.4 keV we derive 
upper-limits of 22 eV and 118 eV on the EW's of the nar- 
row and broad components to the 6.4 keV line respectively, 
which are commonly detected in the spectra of normal 



Seyferts and quasars (Porquet et al. 
let al.||2005| ). 
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Fig. 6. Upper panel: data and best fit model of an absorbed 
power-law to the PN, MOSl and M0S2 (from top to bot- 
tom). The three other panels show the residuals of the fit 
in terms of sigma. 



3.3. Optical data 

3.3.1. XMM-Newton optical/UV monitor 

The optical/UV monitor (OM) detector on board of XMM- 
Newton observed NGC 4278 for six 2880 s exposures in the 
UVM2 and UVWl filters. The central region of NGC 4278 
in the UVM2 filter is detected as a point-like source 
not showing any sign of diffuse emission from the host 
galaxy. The SAS tool omichain produced an OM source 
list based on a 3" aperture and background taken from 
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an aperture with radii 3.7"-4.3". The nuclear corrected 



count rate in the UVM2 is 0.56 counts s 
sponds to a flux density of 1.18 x 10~^^ 



which corre- 



ergs cm 



A- 



(using a count rate to flux conver sion factor of 2 . 103 x 

(|2009l) 



Carduflo et ah 



10 ergs cm ^ A ^ counts ^ ^). 
show the observation of NGC 4278 made with the 
WFPC2/F218 fllter (band-pass: 2000-2500A, comparable 
to UVM2 band-pass: 2050-2450A) in 1994 and in 1995. It 
is clear in their figures that the dominant source of emission 
in this specific band-pass is the central source. Therefore, 
we can consider that the flux calculated from the UVM2 
filter corresponds to the nuclear fiux. We decided not to 
use the UVWl image as it shows strong diffuse emission, 
most likely from the host galaxy, around the central source. 

3.3.2. HST 

NGC 4278 was observed four times with each of the 
ACS/WFC F475W and F850LP fihers between December 

2006 and January 2007. For every filter, we derive the flux 
for two pointings where the source lies away from the CCD 
edge. We use the flux derived from the Hubble Legacy 
Archive (HLA)^° DAOPhot tool, originally in counts s~^, 
and convert it to ergs cm~^ s~^ using the conversion equa- 
tion Fa = C/(e/ JS\dX) where C is the count rate, £/ is 
the fraction of the point source energy encircled within N 
pixels, and JS\dX is the total imaging point source sen- 
sitivity". We find a flux in the F475W fllter of 6.64 x 
10~^^ ergs cm~^ s~^ and 8.02 x 10~^^ ergs cm^^ s^^ for ob- 
servations made in December 23 2006 and January 2 2007 
respectively. The fluxes derived for the F850LP filter are 
7.45 X 10~^^ergs cm~^ s"^ and 8.08 x 10~^^ergs cm~^ s"^ 
for observations made on December 23 2006 and January 2 

2007 respectively. 



4. Spectral energy distribution 

We compile the available photometry data of the nucleus of 
NGC 4278 from the literature and add the data derived in 
this work to construct the NGC 4278 SED. In order to min- 
imize any contamination from the host galaxy, high angu- 
lar resolution is required at all wavelengths. Therefore, only 
data coming from observations with apertures less than 10" 
are used. Table 3 lists the data used to derive the NGC 4278 
SED shown in Fig. [T) 

In radio wavelengths, and in order to minimize the con- 
tamination from the radio jet emission, milli-arcsec data, 
measured through VLBI techniques, are used when avail- 
able. Otherwise, VLA observations with sub-arcsec reso- 
lution are applied. All of the near-IR (1-3 fim) data for 
NGC 4278 should be considered as upper limits as the reso- 
lution of those observations is between 5" and 10". Only the 
10-20 fim mid-IR point could be considered as less affected 
by emission from normal stellar populations and may rep- 
resent to some extent the nuclear emission (thermal emis- 
sion from dust grain can not be ruled out). Optical and 
UV data are obtained from Hubble Space Telescope (HST) 

^ http://web.archive.org/web/20050322143642/ 
xmm.viIspa.esa.es/sas/documentation/watchout/uvflux.shtml 
http:/ /hla. stsci.edu/hlaview. html 

more details can be found at 

http: //www. stsci.edu/hst/acs/ documents/handbooks/cyclel8/| 
c09_expcalc3.html 
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Fig. 7. The spectral energy distribution of NGC 4278. Stars 
derive from non simultaneous radio to UV observations. 
The open circles are the simultaneous XMM-Newton UV 
and X-ray fluxes. The dots represent the contemporary 
HST/ ACS optical and the Chandra (obs. ID: 7081) X-ray 
fluxes. 



observations where 1" apertures are used. One UV data 
point coming from the XMM-Newton OM telescope is used 
(see sect. 3.3.1). X-ray data, measured in this paper, and 
derived from Chandra and XMM-Newton observations are 
used. 

The Galactic extinction t oward NGC 4278 is 0.11 mag 
(de Vaucouleurs et al. 1991). The same value is derived 
using the Galactic hydrogen column density and applying 
the equations E(B-V)=NH/(5.8x lO^i) mag and Ay /E (B- 



V)=3.1. Internal extinction is reported in |Ho et ah ( 1997) to 
be according to a Balmer decrement Ha/H/:(=2!50 lead- 
ing to E(B-V)=0. Infrared, optical and UV data are de- 
reddened according to thes e extinctions and ado pting the 
Galactic extinction curve of Cardelli et al. ( 1989 ) . 



The radio slope between 5 and 22 GHz is somewhat 
inverted (a 0.6), typical for LINERsI/LLAGN (IHo 



1999). Radio loudness is quantifled using the radio to X 
ray ratio deflned as Rx = iyL^{6c m)/Lx, where Lx is 
the luminosity in the 2-10 keV band fTerashima & Wilson 



,2003) . Following this criterion, NGC 4278 is considered as 
a radio-loud LINER with Rx ranging between 2.4 x 10~^ 
and 3.5 x 10"^ where Rx = 3.162 x 10"^ is the bound- 
ary for radio loud systems. There is a tendency for a 
maximum in the SED at mid-IR wavelengths where stel- 
lar population emission is less effective, although ther- 
mal emission from dust grains can not be ruled out. The 
UVM2 flux derived from the OM telescope is used to 
derive an optical to X-ray slope aox and a radio loud- 
ness parameter Ruv where Q;oa:=0.384 log {L2kev/L ^ °) 

and i?,7v=log(L^(6cm)/L^,(2500)A). These two parame- 
ters are later used for comparison purposes with LINER 
and normal Seyfert galaxies. We derive a 2 keV luminosity 
vL^{2 keV) « 4 X 10"'" erg s~^ from the unfolded model 
of the XMM-Newton pn spectrum which gives an aox of 
— 1.2. We emphasize that this value of aox is determined, 
for the first time for NGC 4278, from simultaneous UV and 
X-ray observations. Rjjv is about 2.7 and the 2-10 keV, ax 
(r — 1), spectral index is ^1.2. 
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Table 3. Multiwavelength nucleus data for NGC 4278. 



V 




Aperture 


aatellite/nlter 


Date 


(Hz) 


(ergs s"^) 


(arcsec) 






Lgsxio'** 


7.83x10''" 


3.34 


Ciandra/ACIS-S 


February 2007 


1.93x10^* 


5.56x10*° 


10 


XMM-Newton /EPIC-pn 


May 2004 


1.21x10" 


1.32x10*° 


3.34 


Ciiandra/ACIS-S 


February 2007 


1.21x10" 


5.96x10*° 


10 


XMM-Newton /EPIC-pn 


May 2004 


1.42x10^'^ 


1.05x10*^ 


0.28 


HST/F218W 


January 1995 


1.42x10^'^ 


6.75x10*° 


0.28 


HST/F218W 


June 1994 


1.30x10^^ 


1.21x10*^ 


3 


XMM-Newton /U VM2 


May 2004 


6.32x10^'' 


3.03x10*° 


0.15 


HST/F475W 


January 2007 


6.32x10^'' 


2.51x10*° 


0.15 


HST/F475W 


December 2006 


5.76x10^"' 


1.13x10*° 


1.26 


HST/F555W 


May 1994 


3.65x10^'' 


6.95x10^" 


<6 


UST/F814W 


May 1994 


3.32x10^" 


2.83x10*° 


0.15 


HST/F850LP 


January 2007 


3.32x10^" 


2.61x10*° 


0.15 


HST/F850LP 


December 2006 


2.40x10^" 


4.56x10*^ 


8.00 


1.54 m Mount Lemmon/J 




1.82x10^" 


7.77x10*^ 


7.20 


3.8 m UKRIT/H 




1.36x10^" 


4.50x10*^ 


7.20 


3.8 m UKRIT/K 




8.69x10^^ 


9.57x10*^ 


8.00 


1.54 m Mount Lemmon/L 




2.97x10" 


1.15x10*^ 


5.60 


1.54 m Mount Lemmon/N 




2.20x10^° 


5.36x10^® 


<1 


VLA 


August 2003 


1.50x10^" 


4.00x10^® 


0.15 


VLA 


September 1999 


8.40x10^ 


2.66x10^** 


(3.2xl.9)xl0"^ 


VLB A/ VLA 


August 2000 


5.00x10" 


2.00x10^** 


(3.2xl.9)xl0"^ 


VLBA 


July 1995 


5.00x10" 


2.25x10^** 


2.5x10"^ 


VLB A/ VLA 


August 2000 


1.67x10" 


1.00x10^® 


5.0x10"^ 


VLB! 


May 1979 



Reference 



This work 
This work 
This work 
This work 



Cardullo et al 



Uardullo et al 
'i'his work 
This work 
This work 



( 2009 1 



Lauer et al. 


J2OO5 1 


Capetti et al 

' nil! 


[(2002 



This work 

iHeckman et al. Vl983|) 
Longmorc & Sharpies ( 19821 
Eongmorc & Sharpies ( 19821 
Heckman ct al. (1983 
Heckman et al. (1983 



Giroletti et al., (j2005 |) 



et al 



6iroletti ct al.lWO^] 
Uirolctti ct al. (2005) 
Giroletti ct al. (20C 



[jones et al.,(|198l| 



Note: Infrared, optical and UV luminosities are corrected for galactic extinction. 



5. Discussion 

5.1. X-ray properties 

The NGC 4278 Chandra observations are well fitted with 
a combination of an absorbed power-law and a ther- 
mal component. The temperature of the thermal model is 
0.62 ± 0.04 keV for the different observations. The power- 
law photon index, F, has an average of 2.2lg'2. As for the 
XMM-Newton observation, an absorbed power-law alone is 
sufficient to fit the spectra well with a F of ~ 2.1. The in- 
trinsic hydrogen column density affecting the power-law is 
of the order of 10^" cm~^ in both cases. The fluxes derived 
for the Chandra observations vary between 5.9 x 10""'^^ and 
18.1 x 10~^'^ ergs cm~^ while the flux measured for the 
XMM-Newton observation is 33.8 x 10~^^ ergs cm~^ s~^. 



Terashima & Wilson (2003) studied a sample of LLAGN 
nuclei, that display a compact radio source, observed with 
Chandra. The authors found that the 0.5-8 keV spectra of 
all the sources could be explained with an absorbed power- 
law with most of the objects having a hard X-ray spectrum 
with a photon index <2, including NGC 4278. Our analysis 
of the Chandra and XMM-Newton observations showed a 
somewhat softer 0.5-8 keV spectrum with F greater than 2. 
The NGC 4278 observation treated by [Terashima & Wilson 



(2003) was affected by pile-up at a level greater than 10% 
and at s uch a high level, th e pile-up model of ,DaviS| ( )200lj ) 



used by 
up, is n 



Terashim a fc W ilson (2 0^ to c orrect 
Gonzalez-Martm et al. 



Tor pile- 

(|2009|) studied 



sleading 

a sample of 82 LINERs observed with Chandra and/or 
XMM-Newton and fit the spectra, in most of the 



cases. 



nent. The authors find a median value for the photon in- 
dex, F, of 2.11 ± 0.52 and a thermal component tempera- 
ture of 0.54 ± 0.30 keV, compatible with the values we find 
for NGC 4278. The authors find a bimodal distribution for 
the 2-10 keV corrected luminosities they derived for their 
sample centered at L(2-10 keV)= 10"^^ ergs s~^ and L(2- 
10 keV)= lO-^i ergs s'^. At a lower flux, NGC 4278 lies in 
the low luminosity part of the group with a 2-8 keV lumi- 
nosity of 8 x 10^*^ ergs s"^ while at higher flux, NGC 4278 
moves to the high luminosity group with a 2-8 keV luminos- 



ity of 6 x 10''*' ergs s ^. Additionally, Gonzalez-Martin et al. 



(j2^ 



2009| ) studied one Chandra observation (obs. ID: 7077) of 
NGC 4278 and found a photon index of 2.6^0 3 and a ther- 



mal temperature of 0.5_q ]^ keV comparable to the values we 
derive for the same observation. However, they find an ob- 
served 0.5-10 keV flux of ^4x lO"'^'^ ergs s~^ cm~^ which is 
half the value we derive for the same observation. However, 



Gonzalez-Martm et al. ( 2009 ) did not correct for the ^8% 
pile-up fraction that could explain the discrepancy between 
the two derived flux values as pile-up tempt to underesti- 
mate flux measurements (see Appendix A) . 



5.2. Variability 

At radio wavelengths, NGC 4278 does not show signiflcan t 



inter-year variability at 2 or 3.6 cm (Nagar et al. 2002) 



with a combination of power-law and a thermal compo- 



At 6 cm, [Giroletti et al.| ([2005|) discovered modest vari- 
ability on yearly time-scales with the flux decreasing by a 
factor of 10% b etween 1995 and 2000. At UV wavelengths, 
Cardullo et"aI1 (|2009|) discovered that NGC 4278 is vari- 
able on months timescales showing that its flux increased 
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Fig. 8. The long term X-ray variability of NGC 4278. 



■^2 times between two different HST observations. This is 



( 2005 ) who discovered vari- 



in agreement with |Maoz et al. 

ability in 12 of 16 L1JNEKs1/LLAII]T objects observed with 
the HST ACS/HRC F250W and F330W filters. Nine of 
the 12 galaxies that show UV variability in their sample 
have detected radio cores suggesting a connection between 
UV variability and the presence of radio cores. NGC 4278 
enforces this correspondence as it shows a luminous radio 
core as well as significant UV variability. Additionally, we 
detect optical nuclear variability on years time-scale as the 
fluxes we derive from HST ACS F475W and F850LP filters 
are 2 to 4 times higher than optical fluxes derived from 
HST WFPC2 F555W and F814W filters (we do not dis- 
card optical variability on shorter time-scales but no intra- 
year optical observations were ever made for NGC 4278). 
In X-rays, NGC 4278 show no significant short time-scale 
(hours to days) variability during the different Chandra ob- 
servations with the "excess variance" , cr^^^ , never exceed- 
ing 0.01, with an error < 15%, which means that NGC 4278 
does not follow the same trend as Seyfert galaxi es show- 



ing stronge r variability with decreasing luminosity ( Nandra 
et al.||1997[ ). On the other hand, during the XMM-Newton 



observation, NGC 4278 exhibits short-time scale variability 
with the flux increasing by a factor of 10% on a 1 hour pe- 
riod. Moreover, a significant monthly time-scale variability 
is observed between the different Chandra observations and 
between Chandra and XMM-Newton observations (Fig. Isl) 
where the flux increased by a factor of ^3 on a few months 
time-scale and by a factor of 5 between the faintest and 
the br ightest observation separated by ~3 years. |Ptak et aT] 
1998 ) suggested that the lack of variability on short (< 1 
ay) time-scales and the detection of variability on longer 
time-scales could be due to the fact that the regions re- 
sponsible for the X-ray emission are larger in LINERs and 
LLAGN than in typical Seyfert galaxies. One possibility is 
that the X-ray emission originates from an optically thin 
geometrically thick flow (RIAF) rather than from a geo- 
metrically thin accretion disk. 

This scenario could be present during the different 
Chandra observations corresponding to the lowest flux 
level. The lack of variability on time-scales of hours to days 
found during the Chandra observation could be compatible 
with a large emission region, such as an extended RIAF 
emission. However, in AGN an anti-correlation exists be- 



tween the 2-10 keV X-ray variability and the mass of the 
black hole (Lu & Yu 2001 1 wit h some e xceptions whenever 
the accretion rate is high 



(.e.g., 



Reeves et al. 2002). Hence, 



the lack of variability during the Chandra observations in 
NGC 4278 is consistent with the high black hole mass of 



3 X 10** M0 (Wang & Zhang 



2003 



Chiaberge et al. 2005 1 



and the low accretion rate. 

The short time-scale variability (t~1.5 h), observed 
during the XMM-Newton observation (where the X-ray flux 
level was highest), could be explained by a disturbance in 
the (extended) RIAF or by the presence of a more compact 
region, e.g., an accretion disk truncated at a lower radius. 
In the context of X-ray fluctuations in a RIAF scenario. 



Manmoto et al.| ( 1996[) and [Takeuchi fc Mineshige| ([1997 1 
were able to explain the X-ray s hots seen in Cygrius X- 



light curves ( [Negoro et al. 19941 as due to a disturbance 
in an extendwTTnAF. However, contrarily to Cygnus X-1 
that switches from a soft state to a hard state just after the 
peak, NGC 4278 remains in a soft state during the whole 
observation (including the variation). Additionally, if the 
variation in the case of the XMM-Newton observation is 
the result of a disturbance in the accretion flow, the max- 
imum radius can be determined from the time duration of 



the v ariation from the free-fall time-scale ( Manmoto et al. 
1996|: 



Tiff) « 4 X 



103r„ 



3/2 



M 



lOMp 



(1) 



where M is the BH mass, T{ff) is the free-fall time 
(the duration of the shot), r is the transition radius, and Vg 
is the gravitational radius, GM/c^ . This yields to a transi- 
tion radius of r ~ 1.3 between the RIAF and the classi- 
cal thin disk. Therefore a disturbance in an extended RIAF 
region appears unlikely to explain the short time-scale vari- 
ation observed from NGC 4278, although further works on 
variability from RIAF are required to establish any defini- 
tive conclusion. Therefore, the presence of a more compact 
region during the XMM-Newton observation, (i.e. the ac- 
cretion disk is similar to that observed in low luminosity 
Seyfert 1 galaxies or most probably truncated at a lower ra- 
dius than that during the Chandra observations) appears 
to be a more viable explanation for this short time-scale 
variability. 

5.3. NGC 4278 BED comparison to other LLAGN 

We now turn to examine whether the X-ray continuum 
emission of NGC 4278 originates from a classical optically 
thick geometrically thin accretion disk or from a RIAF 
structure forming around a black hole (BH) whenever the 
radiative efficiency is very low. The bolometric luminos- 
ity Lfto/ is estimated by integrating the SED from radio 
to X-ray data, neglecting upper limits such as points in 
the IR. As mentioned in sect. 5.2, NGC 4278 as most of 
LINERsI/LLAGN, show a high degree of variability on 
months timescales. Hence, two bolometric luminosities were 
derived where the highest X-ray flux [XMM-Newton flux) 
is used to derive the highest Lhoi and the lowest X-ray 
flux (Chandra obs. ID 7080 flux) is used to calculate the 
lowest Lboi- This gives Lbo/ between 2.21 x lO''^ erg s""'^ 
(low state) and 2.88 x 10''^ erg (high state) which leads 
to an Eddington ratio Lboi/LEdd between 5.7x10"^ and 
7.4x10"^. Similar Lboi/LEdd value were derived for the 
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LINER galaxy M81 ( jPian et al.1|2010[ ). Two possibilities 
emerge to explain the faint luminosities in both the low 
and high state; accretion is occurring at sub-Eddington rate 
an d/or the accre tion is radiatively inefficient. 



42 



Maoz 



( 2007] ), when studying a sample of 13 nearby 
galaxies with LINER nuclei, argued that their SED do re- 
semble the SEDs of normal Seyfert galaxies, and hence the 
mode of accretion, "geometrically thin optically thick ac- 
cretion flow", and the X-ray continuum, "the result of in- 
verse compton scattering of UV/soft X-ray in a hot opti- 
cally thin plasma in a corona" , are also similar. The author 
based his arguments on two parameters: the optical to X- 
ray, aox, slope (calculated using the 2500 A flux), and the 
radio loudness, Ruv, (also based on the 2500 A flux). He 
showed that the aox values derived for his sample overlap 
with the values deri ved for a sample of broad-lined AGN 
(Steffen et al. 2006) and demonstrated that even though 
LINERsI/LLAGN are louder by a factor of -100 than high- 
luminosity quasars, most of the LINERs are considered ra- 
dio quiet compared to their 2500 A flux. NGC 4278 
is -1.2 (see sect. 4), comparable to the values derived for 



Maoz (2007) sample of LLAGN and it follows the exten- 
sion to low luminosities of the relation between the aox 
and the 2500 A luminosity derived for a sample of Seyfert 



galaxies ( Steffen et al. 



2006). Moreover, NGC 4278 Ruv is 



2.7, thus it stands in the radio quiet branch of the relation 
between Ruv a nd the 2500 A lum i nosity for a sample of 



luminous AGN (Sikora et al 
that our results may be slight 



2007 Maoz 2007). We note 



y biased by the fact that the 



2500 A luminosity we derive for NGC 4278 comes from the 
UVM2 XMM-Newton optical monitor telescope where nei- 
ther the peak A nor the effective A are as cl ose as t he HS T 



Maoz 



(2007). 



ACS/HRC F250W filter to 2500 A used by 

A RIAF model has been suggested as a solution for the 
low luminosity observe d in a number of LIN ERsI/LLAGN; 
M81 and NGC 4579 dQuataert et al. 



(Gammie et~al] |l999|), NGC 3998 



1999|), 
^PtaTet 



NGC 4258 



al. 



NGC 1097 fNemmen et al.|2006D, and, more rece nt. 



2004|, 

/, for a 

sample of 24 LINER sources (INemmen et al. | 2010[ ). Authors 
of these papers emphasize the fact that SEDs of LLAGN 
in general are in fact different than normal Seyfert SEDs. 
They differ in a number of ways, mainly, the lack of the 
canonical "big blue bump" at UV wavelengths, observed 
in almost all of the luminous galaxies, which is associated 
with thermal emissio n from a geome trically thin, optically 
thick accretion disk plo||2005[ |2002| . Instead, a maximum 
in the SED is peaking m the mid-lR. The presence of radio 
jets is a common feature in LLAGN as for they tend to be 
radio loud systems. Additionally, the lack of a broad Fe Ka 
emission line at 6.4 keV in most of the LLAGN studied up 
until now, originating from a standard disk in typical lumi- 
nous AGN, might suggest the truncation of this structure 
at a radius much larger than the innermost stable circular 
orbit (ISCO) where a radiatively inefficient accretion flow 
(RIAF) structure forms in the inner part. We note that the 
last argument, lack of a broad Fe K line, used to point at 
the existence of a RIAF in LINERs is not decisive since a 
small fraction, 25-30%, of Seyfert galaxies exhibit br oad- 
ened emission in the iron K band ( [Nandra et aI1[2007| . 

The absence of a narrow Fe Ka emission line at 6.4 keV 
with a 22 eV upper limit on its EW is in agreement with the 
possibility of a RIAF presence at the center of NGC 4278. 
However, the large upper limit of 118 eV on the EW of 
a broad line at 6.4 keV is not discriminative. In order to 
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Fig. 9. NGC 4278 SED: black dots represent the radio data. 
Green dots represent the contemporary HST/ AGS and the 
Chandra 7081 observation fluxes. The light -blue dots cor- 
respond to the simultaneous XMM-Newton CM and pn 
fluxes. The average SED of a sample of LINERs is shown a s 
solid black line (figure adopted from Nemmen et al.|[2010 ) . 
The dashed line represents the effect of extinction correc- 
tion. For comparison, the average SED of Elvis et al. (1994) 
for radio-loud (dotted line) and radio-quiet (dot-dashed 
line) quasars. All SEDs are normalized to l^m. 



examine more accurately the possibility of a RIAF struc- 
ture and to identify the accretion mode occurring at the 
nucleus of NGC 427 8 we plotted in Fig. M (adopted from 
Nemmen et al.||2010[ ) fluxes derived for NGC 4278 from si- 



multaneous or close-by observations. The black solid line 
represents the average SED of a sample of LINER nuclei 
(the dashed black line represents the effect of extinction), 
the dotted red line a nd the dot-d a shed b lue line represents 
the average SED of Elvis et al. (1994) sample of radio- 



loud and radio-quiet quasars, respectively. All SEDs are 
normalized to 1/im. Light-blue dots represent the simulta- 
neous XMM-Newton UVM2 flux at UV wavelengths and 
the EPIC-pn fluxes at X-rays. The green dots represent 
the ACS/WFC optical data and the Chandra obs. ID 7081 
X-ray fluxes. The EPIC-pn and the OM observations were 
taken simultaneously and the time difference between the 
Chandra and the HST data is less than two months and 
therefore assumed simultaneous as large optical and X-ray 
variability is observed only on time-scales of several months 
(see sect. 5.2). Black dots represent the radio data and are 
considered constant on time-scale of several years. 

If we compare the NGC 4278 SED during the low 
Chandra X-ray flux (black and green dots combined) to the 
other SEDs shown in Fig. [9j we can consider NGC 4278 as 
a radio loud LINER galaxy. The radio emission is more im- 
portant than the one observed in radio loud quasars, plus, 
optical and X-ray luminosities correspond to the luminosi- 
ties derived for the sample of LINERs. On the other hand, 
the NGC 4278 SED during the XMM-Newton observation, 
where the X-ray flux is the highest, (black and light -blue 
dots combined) seems to be similar to the SEDs of a small 
number of low luminosity Seyfert 1 galaxies, in which th e 



UV to X-ray slope is almost flat (Papadakis et al. 2008) 



where the authors showed that the accretion mode in lo' 



iW 
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luminosity Seyfert 1 galaxies [Lx/LEdd < 10 resemble 
the accretion mode detected in luminous AGN. 

As NGC 4278 shows a high degree of variability on 
months time-scale fluctuating apparently between a low- 
state SED (black and green dots) and a high-state SED 
(black and light-blue dots) where short time-scale variabil- 
ity is detected, we consider that the emission mechanism re- 
sponsible for the bulk of energy from radio to X-rays during 
the low state is jet/RIAF dominated; a thin accretion disk 
most probably truncated but at a lower radius than that 
during the low-state level would be responsible for the op- 
tical to X-ray emission during the high-state SED. 



6. Summary and conclusion 

We have studied in detail seven X-ray observations for 
NGC 4278, six ACIS-S Chandra observations plus one 
XMM-Newton observation. The observations covered a 
three year period allowing an extensive variability study. No 
short time-scale (minutes to hours) variability is detected 
during the Chandra X-ray observations, while during the 
XMM-Newton observation, where the highest flux level is 
observed, the flux increases by a factor of 10% in a 1 hour 
period. A significant months time-scale variability is ob- 
served between the different Chandra observations and be- 
tween Chandra and XMM-Newton observations where the 
flux increased by a factor of ~3 on a few months time-scale 
and by a factor of 5 between the faintest and the brightest 
observations separated by ~3 years. We checked the valid- 
ity of increasing variabilit y with decreasing lu minosity de- 



tected in Seyfert galaxies ( Nandra et al 



1997), NGC 4278, 



like most of LINE R sources, does not follow the same trend 
( Ptak et al.||1998 ). However, based on the anti-correlation 
relation between variability and the mass of the black hole 
with some exc eptions whenever the accretion rate is high 
(Reeves et al.||2002), no short time-scale variation is ex- 



The temperature of the 
Because of the XMM- 



pected in JNGC 4278 as it harbors a ~ 3 x 10® black 
hole and accretes at low rate. However, the short-time scale 
variability detected during the XMM-Newton observation 
might suggest a change in the central engine of NGC 4278 
at that point. No variability is detected in the hardness 
ratio during the observations being soft in all cases. 

The best spectral fit found for the six Chandra observa- 
tions is an absorbed power law plus a thermal component. 
The intrinsic hydrogen column density affecting the power 
law has an upper limit of 6.7 x 10^° cm~^, while the av- 
erage photon index F is 2.2^q2. 

thermal component is 0.63;to'o4 keV. 
Newton moderate angular resolution compared to Chandra, 
we subtracted, in a 10" radius, the contribution from X-ray 
binaries and diffuse emission calculated from the Chandra 
observations. We find that the power-law emission is dom- 
inant over the whole 0.5-8 keV band and therefore no ther- 
mal component is needed. The hydrogen column density is 
(3.8 ± 0.8) X 10^" cm~^ affecting a power law with a pho- 
ton index F of 2.1. No Fe Ka emission line at 6.4 keV is 
detected with a 22 eV upper limit on its equivalent width. 

We measured a UV flux of about 1.2 x 
10-15 ergs cm-2 s^^ A^i from the XMM-Newton 
optical monitor observation of NGC 4278. The advantage 
of such a measurement is the simultaneity with the 
X-ray observation crucial for comparing such sources 
with normal galaxies as NGC 4278 shows significant 



variability on months time-scales. We calculated optical 
fluxes coming from HST ACS observations in two different 
filters (F475W and F850LP). These fluxes are 2 to 4 times 
higher than optical fluxes, coming from the HST WFPC2 
F555W and F814W filters, observed ^ 12 years earlier. 
Our optical fluxes are contemporary to the Chandra fluxes 
coming from the 7081 observation as the two observations 
are separated by less than two months. 

In order to determine the origin of the emission mech- 
anism responsible for the bulk of energy from radio to X- 
rays and to assess the geometry of the central engine in 
NGC 4278 we plotted the SED derived from simultaneous 
observations (radio^^-|-UV-|- X-rays and radio-|-optical-|-X- 
rays) with the SED of a sample of LINER sources 



(Eracleous et al. 2010) along with radio-quiet and radio - 
loud quasars (fig. 8; adopted from [Nem men et al.||2010[). 
As NGC 4278 appears to fluctuate in optical, tJV and X- 
rays on month-timescales, it appears that NGC 4278 SED 
is consistent with LINER SEDs when the X-ray emission 
is weak and therefore we suppose a jet/RIAF origin for the 
radio to X-ray bulk of energy. On the other hand, when 
the X-ray emission is important, the NGC 4278 SED is 
similar to the SED of low luminosity Seyfert galaxies and 
a thin accretion disk is responsible for the bulk of energy 
from UV to X-rays. Hence, we consider NGC 4278 to alter- 
nate between LINER-like and Seyfert-like nuclear activity 
depending on the strength of its X-ray emission. 
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Appendix A: Going through pile-up 

When a source is sufficiently bright that two or more photons are inci- 
dent on the same CCD pixel in a single CCD frame time, the energy 
of these photons will be added and the observation suffers photon 
pile-up. Low energy photons will migrate to high energies creating 
fake high energy events. Furthermore, if the added energy of the pho- 
tons are greater than the on-board spacecraft threshold, the photon 
will be rejected by the spacecraft software. Accordingly, piled-up ob- 
servations will suffer a decrease in the total observed count rate and 
the spectral shape of the source will be distorted. Although there is 
no simple way to overcome pile-up, two methods have proved effec- 
tive. The statistical method of Davis ( 200~l]l reconstruct piled events 
in a spectrum by applying a non-linear integral equation that will 
connect X-ray source spectra to CCD instrument spectra that takes 
into account the possibility of photon pile-up. A second more secure 
approach is to discard spatially and spectrally distorted piled events 
from the cor e of the PSF, and st ay with non-piled events in the wings 
of the PSF | |Broos et al.|[T998| l. This method also reconstruct tem- 
poral features and is applicable to sources with any degree of pile-up 
fractions. 

Pile- up could be estimated using the count rate landing in a 3x3 
pixel island and the CCD readout frame time. A simple converter 
from count rate to pile-up fraction at the on-axis position is written 
by Dr. M. Tsujimoto and is provided as a C code^'^. Unfortunately, we 
found that all of the seven Chandra observations are affected by pile- 
up. To_overcome_yiisjssue, we decided to follow the method prepared 
by [Getman et ak ( |2005 i^^. The idea is to exclude the piled events 
from the center of the FSF using the acis.extract (AE) tool^^. 

We used AE tools to sample the PSF between radii r^^ and rout, 
the latter fixed to 99% of the PSF, ri„ varies between excluding 0% of 
the PSF to excluding 90% of the PSF in steps of 5%, whereas a step 
of 1% is used to sample the PSF between 90% and 95%. Those annuli 
were then converted from PSF fraction to pixels. The idea is to find 
the smallest ri„ that includes the greatest number of counts for later 
spectral analysis without reducing the source flux due to pile-up. 

Since pile-up is related to the surface brightness of the incident 
photon distribution, we created the radial profile of annuli regions us- 
ing the dmextract tool as described in CIAO and calculated the 
surface brightness of every extracted annulus. We discarded the cores 



giving a surface brightness greater than 1.2 X 10 ^counts s ^ pixel 



minimizing the pile-up fraction to 2%. Fig. |A.1| shows the surface 
brightness profile of the observation 707 8 wi th the black horizontal 
line marking the 2% level of pile- up. Fig. |A.2] rcpresents the ffux pro- 
file of the diminishing annuli with the black vertical line marking the 



The 



code, crSpf.c, can 



be 



found 



,http://www. astro. isas.iaxa.jp/~tsujimot/arfcorr. html 



More 



at 



at 



information can Ee found 

http:/ /www. astro. psu.edu/xray/docs/TARA/ae_users_guide/pileup.txt 

http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 
^® Details on how to obtain a radial profile can be found at 
http:/ /cxc. harvard.edu/ciao /threads/radial_profile/index.py.html| 
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Radius of inner circle (pixels) 

Fig. A.l. Surface brightness profile of NGC 4278 during 
the observation 7078. The sohd horizontal line mark the 
surface brightness value that corresponds to a ^2% pile- 
up fraction. The vertical line shows the radius of the inner 
circle that has to be discarded. 



Radius of inner circle (pixels) 



Fig. A. 3. The photon index, F, as a function of ri„. 
Although it is clear that there is a tendency for F to in- 
crease going away from the PSF center, the different values 
are within the error bars. 



0.5 1 1.5 2 2.5 3 3.5 4 

Radius of inner circle (pixels) 

Fig. A. 2. Inferred flux of NGC 4278 during the observation 
7078 as a function of following the technique for a series 
of excluded core circles. The vertical solid line corresponds 
to the ~2% pile-up fraction. The flux increases to reach a 
plateau when the piled up core is discarded. 



2% level of pile-up. We notice an increase in the flux when the pile- 
up fraction decreases and then reaches a plateau after discarding the 
piled-up core. Table A.l gives the parameters of the Chandra obser- 
vations corrected for pile-up. 

Considering the small exposure time of the observation 398 that 
led to a small number of counts, we decided to discard an inner circle 
sustaining 5% pi le-up fraction. This leaves us with 20% of the PSF. 
We show in Fig. |A.3| that the photon index is increasing when we go 
away from the PSF center but the values remain within the error bars. 
We decided not to include this observation when performing detailed 
spectral analysis to avoid any misleading results that could generate 
un-trustful conclusions. 



Appendix B: Surrounding medium around 
NGC 4278 nucleus 

We turn now to study the environment around the nucleus of 
NGC 4278. Two sources, labeled source 1 and source 2 in Fig.[2](here- 
inafter, SI and S2), within a 10"-radius circle around the nucleus were 
detected in the six observations and therefore complete X-ray spectral 



study can be possible. We consider the rest of the medium as diffuse 
emission plus emission from point— like sources detected in five of the 
six observations or less. Source events and spectra were extracted, for 
both SI and S2, from a polygonal shape around the sources centroids 
that encircles 90% of the PSF. Background annular regions wer e con- 
structed in the same manner as for the nucleus (see Sect.jSTTJ. The 
spectrum and the background for the diffuse emission is constructed 
as explained in the AE users guide section 7.1.2. 

We fit the six spectra for each of the two sources, SI and S2, 
simultaneously with a power— law model (Fig. B.l and Fig. B.2, re- 
spectively). Both fit are acceptable with reduced of 1-06 and 
0.8 for 43 and 27 d.o.f for SI and S2 respectively. The fit re- 

" " for SI 



Q+0.7 



suits in a power— law photon index of 1.0 it 0.2 and 1.3_g 4 
and S2 respectively. The 0.5-8 keV absorption corrected fiuxes are 
2.6toi ^ 10-1'' ergs cm'^ s"! and O.etpi x 10"" ergs cm'^ s"! 
for SI and S2 respe ctively. Adopting t he same distance to the nu- 
~ " this results in a corrected 

-1 for SI and S2 re- 



cleus of 16.7 Mpc (To nry et ar]|2001|l , 
luminosities of 8.7 X 10*° and 2 X lO'^* ergs s~ 
spectively. The power-law spectral fit and the luminosities derived 
for these two sources make them compatible with the low mass X- 
ray binary populations in elliptical galaxies (see [Fabbiano 2006 for 
a review). The hydrogen column density of the intrinsic absorption 
affecting the power-law model has an upper limit of ~ 10^^ cm~^ for 
both SI and S2. We note that a disk black-body model can fit the 
spectra as well with an inner disk temperature of 1.2 and 2.9 keV for 
SI and S2 respectively. 

The diffuse emission is fit with a combination of an absorbed 
power-law and a thermal mekal component (Fig. B.3). We find a 
reduced of 1-3 for 136 d.o.f. The power-law photon index, F, is 
1.8 it 0.2 whereas the temperature of the thermal component is 0.4 it 
0.1 keV. The discrepancy of this thermal temperature to the nucleus 
thermal temperature obtained with Chandra inside a radius of ~ 3" 
could be the fact that the spectra derived for the diffuse emission 
includes part of the less-hot interstellar medium at a larger distance 
from the nucleus (10"). The intrinsic absorption derived from the 
diffuse emission fit has an upper limit of ~ 10^^ cm~^. The corrected 
0.5-8 keV fiux is 5.3 X 10"^'* ergs cm~^ s~' which gives a corrected 
luminosity of 17.7 X 10^* ergs s~'. 

The total 0.5-8 keV corrected luminosity of the two point-like 
sources and the diffuse emission is 28.4 X 10''* ergs s~^. This corre- 
sponds to 4% of the nuclear flux during the brightest Chandra ob- 
servation (obs. ID 4741) and to 12% during the faintest one (obs. ID 
7081). Now if we assume that the raise in the flux during the XMM- 
Newton observation is due to a burst in one of the X-ray point-like 
sources, SI or S2, then the flux of these objects should increase, in 
about 8 months period, by a factor of 50 or 200, respectively. 

To check if such an increase in one of the sources is qualitatively 
possible we calculated the Eddington luminosity of SI and S2 assum- 
ing that these sources are neutron star LMXBs with a neutron star 
mass of about 1.4 Mq. This gives a L^dd of about 1.75x10'^* ergs s"-' 
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Table A.l. Pile- up analysis of the Chandra Observations. 



OBS. ID 


Exposure 


Pile-up % 


fin 


'fout 


Net counts 


Extracted PSF fraction 


Corrected count rate 




ks 




Pixels 


Pixels 






Counts s~^ 


398" 


1.4 


11.5 


1.45 


2.42 


87 


0.2 


0.26 


4741 


37 


20 


1.90 


3.88 


1112 


0.1 


0.25 


7077 


110.5 


7.6 


1.60 


4.17 


2446 


0.15 


0.12 


7078 


51 


14.5 


1.98 


3.24 


1272 


0.1 


0.19 


7079 


105 


13.2 


2.00 


3.30 


2477 


0.1 


0.19 


7080 


55 


5.2 


1.43 


3.29 


1202 


0.2 


0.09 


7081 


110 


5.6 


1.43 


3.34 


2679 


0.2 


0.11 



Note: "The value of ri„ for this observation corresponds to a 5% pile-up fraction. 




Energy (keV) 

Fig. B.l. The best fit absorbed power-law model to SI 
spectra. 



which leads to a super Eddington ratio, Lx/^eddi of 250. Therefore, 
we can completely discard the idea of an increase in the X-ray flux 
of SI and S2 that resulted in the high flux observed during the 
XMM-Newton observation. Hence, The raise of the flux in the case 
of the XMM-Newton observation is definitely due to variation from 
the LLAGN X-ray emission. Consequently, SI, S2, and the diffuse 
emission contribute to only 2% of the total EPIC flux that would 
explain the less— needed thermal component in the XMM-JVewton X- 
ray spectral flt as the emission is almost totally dominated by the 
LLAGN. 



Energy (keV) 

Fig. B.2. The best fit absorbed power-law model to S2 
spectra. 
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Fig. B.3. The best fit absorbed power-law plus a thermal 
component to the diffuse emission spectra. 



